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Summary. — Recent progress in approaches to determine the elastic constants 
of solids starting from the microscopic particle interactions is reviewed. On the 
theoretical side, density functional theory approaches are discussed and compared 
to more classical ones using the actual pair potentials. On the experimental side, 
video microscopy has been introduced to measure the elastic constants in colloidal 
solids. For glasses and disordered systems, the theoretical basis is given for this 
novel technique, and some challenges and recent advances are reviewed. 



1. Introduction 

Solids respond elastically to deformations. Hooke formulated the laws of linear elas- 
ticity as early as 1678 in the anagram 'ceiiinosssttuv' which stood for 'ut tensio sic vis'. 
Solid state physics developed the atomistic view of ordered solids, starting in the classical 
realm and leading to quantized lattice vibrations [1]. Non-equilibrium thermodynamics 
focused on the fundamental principles at work in elasticity [2]. Spontaneous breaking 
of a continuous symmetry leads to equilibrium states, which are characterized by finite 
order parameters, and are degenerate in energy. Homogeneous fluctuations, which re- 
store the symmetry and transform one ordered state into another equivalent one, require 
no energy. If these fluctuations are spatially inhomogencous, they cause a small free 
energy increase. It is captured in the elastic free energy which depends on the gra- 
dients of the fluctuations. Elastic constants, or clastic moduli appear as second order 
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derivatives of the free energy. They can be measured in the long wavelength limit of the 
dispersion relations, which contain the eigenvalues of the dynamical matrix as function 
of the wavevector. The elastic moduli are also connected to the correlation functions of 
the symmetry restoring fluctuations (the so-called Nambu-Goldstone modes), which are 
long ranged. In ordered solids, translational and rotational invariance of the underlying 
Hamiltonian function of the particles is broken, and (transverse acoustic) sound waves 
are the Goldstone modes. This general framework is well understood and discussed in 
textbooks like [3], and many important insights are already possible with simple and 
intuitive considerations using mean field theory. 



Colloidal dispersions recently have moved into the foreground for studying elastic 
phenomena. On the one hand, the rheology of colloidal dispersions is a well established 
area, and has provided important information on the macroscopic viscoelastic behavior. 
Rheology, including the recently introduced micro-rheology, where a probe particle is 
strongly forced through a sample, is especially aimed at studying the non-linear response. 
Plastic, an-elastic, and yielding behavior can be investigated [4]. On the other hand, 
microscopy of particle trajectories has recently become possible and provides detailed 
information on the local motion. The size and time scales of colloidal particles are the 
basis for this unique possibility. Thermal fluctuations can be recorded and spatially 
resolved dynamics investigated [5, 6]. Rare events like the ones dominating nucleation 
can be observed. This review is concerned with the theoretical basis for using video and 
confocal microscopy for determining elastic properties in colloidal solids, especially in 
crystals with defects, and in amorphous solids. 



The elastic energy contributions, the corresponding moduli, and already the elastic 
field become complicated concepts, if disorder is present. Topological defects destroy 
the order and are an important topic with well established consequences on material 
properties; they are not discussed in this mini-review, but are covered extensively in 
the textbook [3]. Point defects, like interstitials and vacancies, do not destroy the long- 
ranged order, but lead to an additional hydrodynamic field, the defect density. While its 
appearance and coupling to the other elastic fields, had been well understood from non- 
equilibrium thermodynamics point of view [2], the microscopic formulae for all elastic 
constants in the presence of defect densities had not been known. Recently, they were 
established within density functional theory [7], as will be reviewed here. The effect of 
defect diffusion on dispersion relations will be discussed, if defects are falsely ignored, or 
taken into account properly. As an introduction, a short summary of the more classical 
approach based on the particle pair potential will be given. The second part of this short 
review then concerns the elastic moduli and dispersion relations in colloidal glass, which 
recently could be measured by video microscopy in two dimensions. It is hoped that the 
review will motivate studies also in three dimensions. 
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2. Introduction to elasticity in the context of thermodynamics 

The free energy of a crystal at constant temperature (dT = 0) obeys the Gibbs fun- 
damental form, the first law of thermodynamics connecting infinitesimal energy changes 

(1) dF = -pdV + fidN + h al3 dll a p . 

Besides the familiar mechanical and chemical terms, it includes a term with a stress 
tensor h a p at constant volume V and particle number N times an extensive strain ten- 
sor U a p — Vu a p. The corresponding work done is SW = J h a p8u a pdV ', and u a p = 
5 (V Q u^(r) + V ( gu Q (r)) is the symmetrized strain tensor field. The displacement field 
u(r) describes the coarse-grained displacement of the particles from their unstrained sites; 
its proper definition will be given below. Because we remain close to unstrained equi- 
librium, the linearized form for u a p suffices in the present context. Einstein's sum con- 
vention also is employed throughout the text, that repeated (Greek) indices are summed 
over all spatial directions. 

The elastic constants are defined as second derivatives of the free energy, which shall 
be taken at constant density 



(2) d2F ' V 



dUafidu-fb n du 7 s 



dh c 



Because of this definition, the tensor of fourth rank obeys the Voigt symmetry relations, 
C a pyS — Cp^s — Cafjs-y = CySa/3, which follow as the strain tensor is symmetric, and 
the C are second order derivatives. 

The above free energy is formulated to capture the existence of long-ranged corre- 
lations in the displacement field, which can rigorously be proven using a Bogoliubov 
inequality [8] (see Sect. 3.3 & 4.1), or motivated via the following heuristic argument 
within linear response: The change in the macroscopic distortion Vu in linear order in a 
small applied homogeneous field h a p, is given by the thermodynamic derivative (suscep- 
tibility) \, x — 5(Vu)/5h\h=o, which according to the fluctuation-dissipation theorem is 
given by (brackets (...) denote canonical averaging): 

k B T X = = ^ff^dv 1 (V Q u(r) V>(r')) = ^J<b V 2 (u(r) u(0)>, 

when spatial homogeneity and isotropy are assumed at long distances for the sake of 
simplicity. Continuing in this simplified consideration, and assuming that the remaining 
integral is given by a rapidly varying function, one arrives at 

2/-/-\ ^bT r/ M /__/_\ k B T 1 



V 2 (u(r)u(r')) = -^<5(r-r'), * ( u (r) u(r')) 



where R cx 1/x is some measure of the elastic constant of the material. This heuristic 
argument for long-ranged correlations (for an isotropic solid it could easily be made 
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rigorous solving the equations of linear elasticity theory [9]) is backed- up by a proper 
and general discussion of the correlations of the displacement field. They are given by 
the (inverse) of the dynamical matrix V a p (q) at wavevector q: 

(3) (<5<(q) Su (q)) - VkeTV-^q) and 2? Q7 (q) C c a ^ s q p qs for q -+ 0. 

The small-<7 divergence like 1/q 2 corresponds to the 1/r long-ranged nature in real space. 
Note the subtle difference, that C c in Eq. (3) is taken at fixed defect concentration c, 
while C™ in Eq. (2) is at constant density n; an example of the (sometimes cumbersome) 
considerations required in disordered solids; see Sect. 3.3 and the Appendix. The Eigen- 
values A s (q) (the subscript s denotes polarization) of the dynamical matrix 2?(q) can be 
considered the spring constants of the system, and are called 'dispersion relations'. In 
atomistic solids, they determine the eigen-frequencies of lattice vibrations, but in col- 
loidal solids the motion is over-damped due to the surrounding solvent. Nevertheless, 
the dispersion relations characterize the elastic energy stored in plane wave distortions 
(with wavelength 2ii/q) of the solid. 

3. Crystalline solids 

This section focuses on microscopic expressions for the tensor of elastic coefficients 
and the dispersion relations in real, mono-domain crystals, pointing to the corrections 
from local disorder or non-affinc deformations. Because of the advances in microscopy of 
colloidal solids, the non-ideality of samples obtained by self-assembly plays an important 
role in the analysis of data and, also, can be studied in detail. Theoretical concepts 
incorporating disorder will be presented, even though, for lack of space, specific results 
on defect diffusion and defect-strain coupling will not be discussed. Those aspects can 
be found in the literature which will be mentioned. 

3T. Born term. - The reference result for the calculation of elastic constants from 
the particle interactions was obtained by Born [1] , and follows from an expansion of the 
potential energy for small displacements of the particle positions around their equilibrium 
sites (i — 1, . . . , N counts through the particles) 

(4) C® M5 = ^ J2 *crr(h3)Rij,f)Rij,6 , 

with <f> Q7 (i, j) = Vq ^V 7 2 ^<l>(Rij), and the actual pair potential $. The derivatives are 
evaluated at the equilibrium positions Ri, and Hij — Hi — Hj. While in ideal crystals, 
which is the situation originally considered, periodic lattice sites were used to define 
the Hi, in later application to colloidal solids the time averaged centers of the particle 
trajectories were employed [5, 6]. Regions without point defects were studied, because 
defect diffusion would lead to large excursions of single particles from their lattice sites, 
invalidating the expansion; see Sect. 3.4. Also, they would destroy the regularity of 
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the real space lattice, so that the dispersion relations would not be periodic from one 
Brillouin zone to another. 

3'2. Fluctuation term. - The above Born term neglects thermal fluctuations, as only 
the change in potential energy is considered. At finite temperatures the free energy needs 
to be considered and its variation with strain. Motivated by the need to analyze their 
atomistic computer simulations, Squire, Holt, and Hoover derived expressions for the 
isothermal elastic constants valid at finite T and for standard particle pair potentials 
[10]. The relations between different expressions were further clarified [11], and the 
formula combining Born and the so-called fluctuation term nowadays are the basis for 
studies of elasticity in solids, including in disordered ones, using computer simulations 
[12]. In principle, they could also be evaluated by microscopy in colloidal solids if the 
pair potential is smooth and well known. 




Fig. 1. - The left two panels show an affine shear deformation of a unit cell of a crystal, where 
all particles move according to the macroscopically enforced distortion. The right two panels 
include an interstitial particle (or the second particle in a basis), whose displacement after the 
affine deformation remains open. Non-affine local particle displacements are especially important 
in crystals with complicated unit cells or in disordered solids. 



Figure 1 presents the conceptual problem when determining elastic energies from 
macroscopically imposed deformations, especially in crystals with complicated unit cells 
or in disordered systems. The imposed macroscopic deformation does not affincly deter- 
mine the individual particle displacement as assumed in the derivation of the Born term. 
Particles can be displaced non-affinely, which even at vanishing temperature causes a 
correction to the Born-term in the elastic constants [11, 12]. The nature of the displace- 
ment of the particles needs to be determined and enters as constraint on the variation 
when calculating elastic constants. In a deformation, the particles move to positions 
where the forces cancel. When considered at vanishing temperature, this condition re- 
quires a 'relaxation' of the particles to the adjacent (potential) energy minima. In the 
canonical ensemble, where the average energy is given, the local stresses also perform no 
work on average, which states that the constraint of vanishing forces also holds at finite 
temperatures. Consider the free energy F — —ksThxZ, which is given by the partition 
sum of the canonical ensemble Z — J dp N dr N e~^ H , taking as Hamiltonian function 

2 

H = Si Im + I Si^j ®( r ij)- The stress field t is the field thermodynamically conjugate 
to strain 

_ 1 OF _ -k B T dZ 

(5) taf> M ~ Vdu^ VZ du aP 
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For performing the derivative, an affinc mapping of arbitrary points in the material 
is applied to the Hamiltonian function, the linear change in Z obtained (note that V 
is varied in this procedure), and the mapping is inverted to obtain the corresponding 
(macroscopic) strain u a p entering thermodynamics. This result for the average (Cauchy) 
stress can be obtained from averaging the Irving-Kirkwood expression of the microscopic 
stress tensor 



(6) i a p = -y 



9 2.^ 



f-f m 2 ^r 1 . \ ru dn 

i=l 1^2 



i ij 



The first contribution comes from the kinetic energy, the second from the Virial and the 
potential energy. The microscopic stress tensor t also determines the conservation law, 
dtja + V ' f)i a p — of the momentum density j(r, t). In unstrained thermal equilibrium, 
the stress tensor averages to the negative of the pressure (see the Appendix): 

(7) t afj \ u ^ s=Q = -p 5 af3 unstrained. 

The elastic constants at vanishing external strain follow as second order derivatives, 
repeating the above procedure, and contain three terms 

dtg/j d {If , 3N . 3NT C -0H\ 

Ca ^ 5 ~ d^s ~ d^ 5 {z J dp dq ta(ie ) 

(8a) = 2nkBT(8 a ~ l 5ps + S a g5p^) 

(8b) + (C^ 5 ) 

(8c) - (<*a/»*7«> - (*a/3><*7«>) 

Term (8a) is kinetic in origin, while the canonically averaged Born term (8b) follows from 
the derivative acting on t. The variance of the Irving-Kirkwood stress tensor, the first 
term in (8c) follows from the variation of the Hamiltonian function in the average, while 
the square of (ty$) originates from the normalization l/Z. In the Appendix it is shown, 
that Eq. (8) gives the elastic constants at fixed, possibly vanishing, defect concentration. 
Because no mapping of the particle positions to periodic lattice sites is involved, the 
above expression has turned out useful for the study of glasses and other disordered 
materials especially at low temperatures; for a more in depth discussion and results on 
glasses see the review [12]. 

3'3. Density functional approach. - Density functional theory (DFT) is especially 
powerful in describing inhomogcneous states, like crystalline solids. Thus there is a long 
body of work, using DFT to obtain the elastic constants of periodically ordered solids, 
which partially is reviewed in [13, 14, 7]. Yet, application of DFT to this problem requires 
an adhoc ansatz of the inhomogeneous density field as expressed by the coarse grained 
fields of elasticity theory Different choices have been made in order to incorporate de- 
fect density, because no fundamental criterion had been known to motivate the ansatz a 
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priori. Also the ansatz was extended to glasses [14], which appears questionable in view 
of the material reported in Sect. 4.3. The recent developments on elasticity in non-ideal 
cystalline solids overcame this problem by going back to the microscopic inhomogeneous 
density field and studying its time dependence [7]. As will be discussed in Sect. 3.4, at 
present the most stringent connection between different approaches to clastic coefficients 
appears to be via the time-dependent equations of motion. Here, before entering into a 
discussion of these finer points concerning DFT, the elastic free energy shall be derived 
from DFT following the standard steps in order to familiarize the reader with the ap- 
proach [15]. DFT expressions are especially suited to hard spheres, which are the basis 
of most colloid-particle models, because the explicit potential and its derivative (in the 
Virial of Eq. 8) are circumvented; these would not be well defined. 

To obtain the second order derivatives of the free energy with respect to its thermo- 
dynamic variables, one starts from the second order change in free energy AJ" due to 
a variation in the density distribution Sp(r) around the periodic crystalline equilibrium 
density [16]: 



where c(ri,r 2 ) is the direct correlation function of a periodic crystal. The average 
periodic density is n(r), which can be observed in e.g. Bragg scattering. The Sp(r) is 
the variation in the microscopic density profile, which varies on length scales extending 
down to e.g. particle size and lattice constant. Assuming that the integral is dominated 
by a smooth variation described with the ansatz 



one can simplify Eq. (9) in order to evaluate the thermodynamic derivatives appearing in 
Eq. (2). Here u(r) is the coarse-grained displacement field, whose symmetrized gradient 
is the macroscopic strain field. The ansatz for the density change is motivated by a 
macroscopic consideration where the smoothly strained average density profile n(r + 
u(r)) is expanded in linear order in u. In the ansatz (10), the fields live on two very 
different length scales. The coarse-grained displacement field describes long wavelength 
fluctuations arising in elasticity theory. The local density profile varies rapidly from one 
lattice plane to the next, as docs the direct correlation function, which is the second 
functional derivative of the (interaction) free energy functional in DFT. Both locally 
varying functions are connected by an important relation in DFT, which is the crucial 
ingredient in deriving the proper thermodynamic derivatives from the above AJ 7 . The 
integrals can be manipulated with the help of the Lovett, Mou, Buff, and Wertheim 
(LMBW) equation. 



(9) 




(10) 



Sp(r) 



u(r) • Vn(r) 



(11) 
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It is a statement of the translational invariance of the underlying Hamiltonian function of 
the system, which is spontaneously broken in the crystal state, but constrains the direct 
correlation function and the average density profile [16, 17, 18]. To employ Eq. (11), 
one enters Eq. (10) into AJ 7 , and expands one slowly varying displacement field ^(r^) 
around ri, the position of the other. The zero order term of AJ 7 vanishes, because of 
Eq. (11) and the first order term does not contribute as c(ri,r 2 ) = c(r 2 ,ri). Stopping 
the expansion at the squared gradient level, because the hydrodynamic variable u(r) is 
assumed slowly varying, one obtains an expression which is quadratic in Vu(r) 



A.F = 



2 



(12) 



r 1 1 d ^ lddr2 (~!^J ~ c(ri,r 2 ))u Q (ri)u lS (r 2 )V a n(ri)V / 3n(r 2 ) 
f f 

■ / d < Vid f V 2 V Q n(ri)c(ri,r 2 )V (3 n(r 2 )u Q (ri) 
(up(r 1 )-up(r 1 )+ V 7 u ( g(r 1 )r 12 , 7 - ^V 7 V <5 u ) g(r 1 )r 1 2 i7 r 1 2 i( 5 + . . . ) 
i J d'Sr A Q ^ 7 , 5 V 7 u Q (r)V < 5U / 3(r), 



with expansion coefficient tensor 
k T f f 

(13) A Q( 3 7(5 = -^y- / / d d r 1 d d r 2 (V«n( r i)) c(ri,r 2 ) (V (3 n(r 2 )) r 12n r 12 ,s 

To derive the required Voigt symmetry of the elastic constants, requires the rotational 
analog of the LMBW equation as explaind in [7]. It arises from the isotropy of the 
underlying Hamiltonian, and leads to (in d — 3 dimensions) 21 independent A-clcments, 
which obey: \af3~ ( s = ^/3a~ ( 5 = ^a/3S~{ — A 7< 5 Q( 3. Then the gradients of the displacement 
field can be symmetrized to get the strain fields, and the four equivalent versions of 
Eq. (12) can be combined. Exploiting the symmetries of A further, the expressions for 
the elastic coefficients can be condensed [19]. Finally, the derivatives in Eq. (2) can be 
evaluated for homogeneous strain to find: 



(14) — \ al j3& + A ( g 7Ct 5 — A Q( g 7 ,5 . 

The C n inherit the Voigt symmetry from the A even though different index combinations 
appear. To evaluate the DFT result for the elastic constants, requires an approximation 
of the direct correlation function of the crystal and for the average density profile. Both 
could e.g. be taken from Ramakrishnan-Yussouff approach to hard sphere crystals based 
on the Percus-Yevick approximation [3]. 

The above result holds neglecting the coupling between strain and density fluctu- 
ations, and had been generalized to include defect densities in various ad hoc ways 
[13, 14, 7]. In order to find the proper way, it is advantageous to combine DFT with the 
Zwanzig-Mori (ZM) equations of motion of the microscopic density fluctuations [20, 3]. 
Thereto, only the reversible couplings are required, which can be formulated without 
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approximations. The criterion for the selection of the variables in the ZM approach is 
given by Bogoliubov's-inequality, see Eq. (20), which yields all long-ranged correlations 
in a crystal [8]. This enables one to generalize the ansatz (10) to include the defect 
density, Sc(r) [7]: 

(15) Sp(r) = -V • (n(r)u(r)) - ^<5c(r) , 

Again, the coarse grained fields, u(r) and Sc(r), are assumed to vary much more smoothly 
than the average crystalline density profile n(r). The first term, which is the only one 
present in an ideal crystal, now properly accounts for the strain-density coupling: the 
shape fluctuation of Eq. (10) is complemented by a volume term, so that the density 
change is given by a divergence-term [13]; this follows from mass conservation as lattice 
motion is identical to particle motion in an ideal crystal. The second term in Eq. (15) then 
couples in the independent defect motion. The ansatz (15) that the slow contributions 
in the microscopic density fluctuations Sp(r,t) are dominated by four coarse-grained 
fields Sc(r) and u(r) solves the (reversible, isothermal) ZM equations for the microscopic 
density fluctuations Sp(r,t) on all length scales. While its uniqueness could not yet 
be shown, it appears unlikely that the DFT approach could lead to another relation. 
Another advantage of the approach is the possibility to invert the relation and obtain 
microscopic definitions of the fields of elasticity theory from the microscopic density 
fields; see Sect. 3.4. The above coarse-graining can be repeated including all elastic fields 
of the free energy density f{u a p,c) = F/V. The / is the starting point for deriving 
hydrodynamic equations in the framework of non-equilibrium thermodynamics [2], and 
(for completeness) is derived from the free energy Eq. (1) in the Appendix. Its term Af 
quadratic in the thermodynamic fields contains all constants of elasticity of a general 
crystal, v, Pap, and C a p-ySi in the general case; see [15] for definitions and calculation: 

A t 1 / V X 2 + x 

A / = o^ -2^ + 2 -5cu a p 

(16) + [C2p y g + v6 a p8 yS + fJ'ap6-y5 + fiaPlJ»y6] u a/3 Uj5 j. 

Because of the coupling between density and the trace of the strain field, there is an effect 
on the elastic constants. The angular bracket [. . .] contains C c , the elastic coefficients at 
constant defect density, while C" is at constant density. 

3'4. On the comparison of the different approaches. - In ideal solids, where each site in 
a periodic lattice is occupied precisely by one particle, defining displacement and strain 
fields and the corresponding elastic constants is straightforward. With point defects, 
like vacancies or interstitials, defining the displacement field becomes more complicated; 
see Fig. 2. Also, the non-affinity of the local particle displacements becomes more im- 
portant [12]. While in equilibrium hard sphere crystals, the defect density of vacancies 
/ interstitials is low at melting (c v = 10~ 4 and q = 10~ 8 , respectively), it increases 
dramatically when the particles have a polydisperse size distribution as is always present 
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Fig. 2. - Schematic two dimensional drawing of a crystal containing point defects, in (a) an 
interstitial and in (b) a vacancy; from [7]. In (c), the definition of the displacement vector Ui of 
particle i is given as difference between actual particle position r< and reference site R^. In case 
of point defects, particles can not be labeled by the lattice positions, as implied when evaluating 
the Born term (4) as sum over lattice sites. 



in colloidal solutions; then the density of interstitials reaches up to Cj = 0.02, before the 
crystal melts for even higher polydispersity [21]. Then, the coupling of defect density 
to strain fluctuations can not be ignored anymore, and the terms /j, and v in Eq. (16) 
matter. Thus it is of interest to compare the different expressions for C a p~i5, to ascertain 
in which situation the correct result is obtained. 

A first observation concerns the Born term Eq. (4), which can be recovered from the 
DFT result Eq. (13) using the so-called random phase approximation (RPA) c(ri,r 2 ) = 
— /3<j>(ri, r 2 ) and appropriately coarse-graining. Reassuringly, this most simple DFT 
approximation to replace the direct correlation function with the true pair potential 
[16], recovers the approximation following from affinc particle displacements at vanishing 
temperature. 

The wave equation, and the complete equations of motion of elasticity theory, provide 
at present the crucial test for the derived relations of the coarse-grained fields and their 
elastic coupling coefficients. For simplification, considering only the wave equation for 
the momentum field, this requires to establish 

(17) mn dt5j a (q,t) = -V a p(q)5jp(q,t) , 

with the dynamical matrix 2? Q| g(q) defined in Eq. (3) from the long-ranged displacement 
fluctuations via the equipartition theorem. The momentum density j(q, t) is the familiar 
one, also entering hydrodynamics. From the phenomenological approach, the wave equa- 
tion and actually the complete long-wavelength coupling including defect density is well 
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established and contains a number of subtle relations like the Voigt symmetries [2, 22]. 
Any microscopic approach needs to recover this. In the DFT approach, the ZM equa- 
tions show that C c and not C™ determines the sound velocities, viz. the long-wavelength 
limit of the wave equation [7]. Beyond this, the dispersion relations for finite q can be 
compared. This route to establish microscopic expressions of elastic constants appearing 
in states with spontaneous symmetry breaking is reminiscent of the comparison of the 
Triczcnberg-Zwanzig and the Kirkwood-Buff expression of the surface tension. These 
two relations, one containing the direct correlation function and the other the actual 
pair potential, can also be identified via the dynamic equations of motion [16]. Reassur- 
ingly in the present consideration, the DFT approach leads to the complete reversible 
and isothermal equations of crystal elasticity with all reversible couplings and their sym- 
metries derived [22]; the generalization to include heat transport and dissipation are 
presently investigated. It would be desirable to establish rigorously that the expression 
(8) including the fluctuation term agrees with the DFT expression (16) containing the 
defect corrections. But while I expect this to hold, the coupling of defects has not been 
addressed explicitly in the former approach. 

To give an example of the comparison of different approaches, Fig. 3 shows the effects 
of local defects on the dispersion relations obtained in computer simulations [23] . The 
eigenvalues of V a p (q) as functions of wavevector, are presented for a hard sphere crystal 
with only few defects. To be precise, Monte Carlo (MC) simulations of 32003 'hard 
sphere Yukawa particles' were performed mimicking the colloidal system of Ref. [6], and 
4 vacancies and 7 intcrstitials were introduced by hand. Equilibrating the defects for 
7 10 4 MC steps, and then measuring the displacement fluctuations over a 'time span' 
At of 100 or 500 MC steps, the dispersion relations were obtained from Eq. (3). Figure 
3 shows that defect motion affects the measurement of the elastic relations depending 
on the definition of displacement fluctuations used. Even a tiny number of defects, 
can destroy the accuracy of the dispersion relations if the defects are mobile and the 
measurement averages over a finite time; this happens in the upper row of Fig. 3. It 
holds in the method defining particle displacements as difference vectors of the particles 
momentary position from the time-averaged center of the trajectory; see Fig. 2(c), and 
Eq. (22) below for the precise definition. This method thus neither requires a periodic 
lattice, nor fails in assuming an one-to-one mapping of particles to (average) sites. Yet, 
it is incorrect and fails because extending the averaging time span At, particles hop on 
vacancy positions and move arbitrarily far following the diffusion of the defects. A (low 
but) finite density of particles appears to be displaced far wider than corresponding to the 
correct local elasticity. In the long run, the variance of displacements of a low but finite 
density (closely related to the defect density) of particles will become arbitrarily large, 
and thus the apparent dispersion relations will be incorrectly low. Note the decrease in 
the so measured dispersion relations in the upper right side panel of Fig. 3. Dashed lines 
are T = calculations based on the potential evaluated at the lattice sites R^. Because 
of the periodicity of the , the theoretically calculated dispersion relations are periodic 
in reciprocal space; note that this does not hold in the simulated ones, because the time- 
averaged centers of the particle trajectories do not lie on a perfect lattice. This, however, 
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Fig. 3. - Dispersion relations of a defective colloidal hard sphere like crystal (of fee structure) for 
various high symmetry directions q in reciprocal space (q given below the x-axis, starting from 
q = 0, the T-point, and symmetrically extended into the 2nd Brillouin zone); colors give different 
polarizations. The results are obtained from MC simulations of the variance of displacements 
according to Eq. (3) [23]. The upper row uses time averaged lattice sites to calculate particle 
displacements, the lower row uses the DFT-expression for the displacement field obtained from 
density fluctuations. Left panels average displacements over 100 MC steps, right panels over 
500. While the DFT dispersion relations improve with better statistical averaging, the upper 
row results neglect defect motion and deteriorate with averaging time. Dashed lines give T = 
calculations using the pair potential based on Eq. (4) evaluated at periodic lattice sites. 



should not be considered an error of the simulation method, because the one-to-one 
mapping of particles to lattice positions Ri is impossible in the presence of point defects. 
The perfect lattice and the periodicity of the dispersion relations are an idealization only 
possible in an ideal crystal. 

The lower panels in Fig. 3, were obtained with displacement fluctuations calculated 
via the inverted relation (15) in the DFT approach. It can best be inverted in Fourier- 
space and expresses the coarse-grained displacement in terms of microscopic density 
fluctuations, e.g. obtainable from a time-dependent DFT approach, [7, 24]: 



(18) 



g 



Elastic properties of colloidal solids with disorder 



13 



with a normalization matrix Af = |?ig| 2 gg- Only spatially modulated density fluctu- 
ations 5p g (q, t) = Jd d r e~*( g+q )' r p(r,t) — n g VSqfi need to be measured, which require 
no definition of a lattice in real space. The crystal structure enters in Fourier-space 
via the reciprocal lattice vectors g and the Bragg-scattering amplitudes n g ; the latter 
are the order-parameters of the crystal, and describe the modulated average density 
n(r) = J2 S e~ 4g ' r n s . Local defects lower the magnitude of the n g (increase the Debye- 
Waller factor, like thermal fluctuations also do) , but do not destroy the reciprocal lattice 
[3]. Because this so-defined displacement is a continuous field (depends continuously 
on r), its dispersion relations are not periodic in reciprocal space; this can be seen in 
Fig. 3, and should not be taken as a defect of the definition. Averaging longer over time, 
increases the accuracy in the sampling of the variance of displacement fluctuations and 
thus reduces the noise in the curves; this can be seen from comparing the lower right 
panel with the lower left one. The temperature in the simulation is quite low, so that the 
T = O-calculation describes the dispersion relations quite well in the first Brillouin-zone. 
Temperature dependence can be seen only close to the Brillouin zone boundary, e.g. at 
and beyond point K, which lies around 1/3 of the T-r distance along the [1,1,0] direction 
[6]. In Figs. 3 and 4, a is the average particle distance. For application of the sums in 
Eq. (18) to real experiments, it is good news that only 14 Bragg peaks were included in 
the results of Fig. 3. 

Reassuringly, the general definition of the displacement field in Eq. (18) reduces to the 
expected low temperature result in ideal crystals. Let particle i fluctuate around its site 
Hi in a periodic lattice. Then up to linear order in the displacement Ui(t) — r^(t) — 
the density fluctuation becomes 

N N 

5p g (q,t) = £ e-'ts+^W » ]T e'^ (l - ig ■ u,(t) + 0((g • u,) 2 )) , 

i=l i=l 

where e _lg ' Ri — 1, and 5 « j was used. Inserting this into Eq. (18), and using that 
J2 S n s g = by inversion-symmetry, gives 

N 

(19) u id - c ^ st -(q,t) = £ e"^ (ui(t) +0((g • u,) 2 )) , 

»=i 

where the normalization by TV could be canceled. This is the expected result for an ideal 
periodic crystal. It is periodic in reciprocal space, u ld cryst - (q + g, t) = u lcLcryst - (q, t), 
which the original expression (18) was not; the difference is hidden in the neglected 
terms of order C((g • u^) 2 ). 

4. Glasses 

This section focuses on microscopic expressions for the tensor of elastic coefficients 
and the dispersion relations in amorphous solids, especially in colloidal glass. Video 
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microscopy measurements of a binary glass former in two dimensions will be presented 
and the theoretical understanding reviewed. Again the question on how to define a 
displacement field and to use the equipartition theorem Eq. (3), will be central. 

4'f. Modified density functional approach. - The first theory to be considered, is 
Masters' DFT-based approach [14], where the coarse-graining of the DFT free energy 
described in Sect. 3.3 is applied to glass. Reconsidering the developments leading from 
Eq. (9) to Eq. (13), all steps appear to go through in amorphous solids as well as in 
crystalline ones; also the equipartition theorem Eq. (3) seems valid at first inspection in 
glass. Yet, the expression (10) for the coarse-grained density field implies a long-ranged 
pair correlation function, and thus a small-g divergence of the static structure factor in 
glass [14]. This should come as no surprise, as Wagner showed that the structure factor 
in crystals obeys the Bogoliubov inequality [8, 7]: 

(20) iWe n)|-)> (e+ ^r i%l ' F ^-'- 

where R abbreviates one of the elastic coefficients depending on g-direction. This 'dif- 
fuse scattering-background' diverges for wavevector q — > 0, but the divergence has a 
non-vanishing amplitude only around a finite reciprocal lattice vector, g ^ 0, where the 
Bragg-peak resides. (The latter arises from constructive interference and thus is an order 
in particle number TV larger than the diffuse background.) The divergence is strongly re- 
lated for all Bragg-peaks, at it is coupled to the long-ranged displacement correlations in 
Eq. (3). In the amorphous glass, the concept of an inhomogeneous average density field 
n(r), whose gradient has to be non-vanishing in the ansatz (10), and which also appears in 
the long-ranged pair correlation function, (g(ri, r 2 ) — 1) — > (Vin(ri)) |r°" 8 r 2 | (^2^(1*2)), is 
quite unclear. First, how should the amorphous average density (a 'one point function') 
be inhomogeneous? Second, why are no long-ranged correlations in static two-point 
functions observed? It seems not easily possible to apply the concepts of DFT and 
equilibrium inhomogeneous density fields to glass, which is a metastable state. Scatter- 
ing experiments should detect the long-ranged density correlations accompanying this 
interpretation of glass as symmetry broken state [14], in the same way as the diffuse 
background diverges close to Bragg peaks in crystals [25] . 

4'2. Replica theory approach. - Another suggestion on how to derive the elasticity of 
glass from an interpretation of glass as a symmetry broken state was made by Szamel 
and Flenner [26]. They showed that in the mean-field replica approach, long-ranged 
displacement correlations, viz. Eq. (3), and (shear) rigidity, viz. finite elastic constants 
(2), emerge from breaking of replica symmetry. In the replicated system, a glass shows 
up in non-trivial two-point off-diagonal densities rijj(r), for different replicas i 7^ j of 
the system. This overcomes the first problem in Sect. 4.1, and keeps the average density 
n(r) = n constant in glass. The long-ranged density correlations, expected follow- 
ing Wagner's Eq. (20), then appear in four-point replica off-diagonal densities. This 
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eliminates the second problem in Sect. 4.1, and predicts finite correlations of the mea- 
surable densities. Of course, glass is an isotropic solid, so that the tensor of elastic 
constants simplifies to the expression containing the two Lame coefficients A and /i only, 
Cap-yS = \6 a p8~,s + (^(Sa-yS p$ + 8a5$f}-y) j where [i vanishes in a fluid. Szamel and Flenner's 
final equation for the shear modulus /j, closely resembles Eq. (13), yet replica off-diagonal 
densities, m^j(r), and direct correlation functions, c^jk^l ( r i , ■ • • , 1*4) , enter, which can 
be evaluated with the replica Ornstein-Zernike equation [26]: 

M = Jdr 1 ...dr i y 13 — — — — [Ci io(ri2 ; r 3 4) - Ci02o(ri2, r 3 4)J • 

This approach thus achieves to combine the concept of symmetry breaking, spatially 
correlated Goldstone modes, and elastic constants with the observed absence of these 
phenomena in the simple one and two point static structure functions of glass. The con- 
nection to the replica theory for the formation of a polymer gel [27] appears interesting, 
as there the order parameter also is a second moment of a static density fluctuation which 
itself has vanishing average. Comparing both approaches might provide additional infor- 
mation on the connection between Goldstone-modes and elasticity in amorphous solids 

4'3. Mode coupling theory approach. - A practical approach to the moduli and elastic 
dispersion relations in glass is inspired by the one for crystals with immobile defects. To 
determine the displacement field using the center of mass of the particle trajectory as 
lattice site [5]; see Fig. 3 upper row for the crystal application (with mobile defects). The 
trajectory has to be averaged for a sufficient long time window but still short compared to 
the structural relaxation time r Q in order to test glassy not fluid- like behavior [28]. The 
theoretical analysis supporting this approach, derives a quasi-equilibrium description of 
the non-ergodic glass state, which is based on the assumption of a kinetic (not thermody- 
namic) glass transition. In this sense, the approach is inspired by mode coupling theory 
(MCT) [29], even though no explicit mode coupling approximations are involved. 

The defining relation which introduces the collective displacement field u q (i) is 

(21) u q (t) = v q (t) . 

Here, the velocity field v q (t) = (1/y/N) X^ili e iq ' r< (*) r,(i) is well familiar from fluid 
hydrodynamics. (Note the convention to normalize fluctuations by 1/y/N, which is used 
in Sect. 4.3. The relation to the momentum density of Sect. 3.2 is j = y/Nmv.) Central 
question is how to integrate Eq. (21) in time. Let time averaging be denoted by an 
overbar so that t\ is the average position of particle i G [1,-^V] during the time interval 
At. Then a possible ansatz for the particle's displacement is Uj(i) = r^(t) — 1\, and 
Eq. (21) leads to 
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Fig. 4. - The 'dispersion relation' of 2D binary glass of soft repulsive colloids for different rescaled 
inverse temperatures F as extracted from the displacements of particles from their time-averaged 
positions [28]; averaging time At w 10 4 s <C r a short enough, so that starting from slightly below 
the glass transition T c « 200, elastic behavior can be seen at intermediate times. (The average 
particle distance is a w 21/xm.). Filled and empty symbols represent longitudinal and transverse 
waves, respectively. The right panel shows typical particle trajectories of a big (green) and a 
small (red) particle at F = 242, shifted to center around the origin. 



to order O(q-Uj). Yet, does this result make sense? Do the displacements remain 
bounded, so that Eq. (22) remains valid for a sufficient window in time? Goldhirsch and 
Goldenfeld solved Eq. (21), differently and used u,-(t) = Yi(t) — rj(0), where time t = 
is some arbitrarily chosen initial time [30]. They considered granular systems and were 
not concerned that the equilibrium-like equipartition theorem (3) does not follow from 
their approach; rather (u* u q ) = (u q (f = 0)* u q (t = 0)) = follows. Yet, the discussion 
in Sect. 3.4 showed that Eq. (22) fails in non-ideal crystals because of defect diffusion; 
recall that it was used in the upper row of see Fig. 3, and gave results deteriorating with 
At, the time span over which the trajectory was averaged, and over which 0(q-Uj(£)) is 
required to remain small. 

Under the assumption of an ideal glass transition (as described by MCT), it can 
be shown that the variance of displacements remains bounded using the Zwanzig-Mori 
equations of motion (EOM). This justifies Eq. (22) as leading order approximation to the 
solution of Eq. (21). An approach, considering glass as non-ergodic state, best starts in 
the ergodic state, and considers simple and well defined correlation functions, which can 
become non-ergodic, viz. take finite values at infinite time. The following mean squared 
collective displacement function is a straightforward collective generalization of the single 
particle mean squared displacement: 



(23) C(q, f) = (Au;(t) Au q (i)) = 2 fdt' (t-t') (v q (t') v q (0)) , 

Jo 
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where the displacement difference is Au q (i) = J dt' v q (i'). Its overdamped EOM 
appropriate for colloidal particles neglecting hydrodynamic interactions can be found 
from standard Zwanzig-Mori steps [28]: 

(24) C(q, t) + [ dt' G(q, t - t') C(q, t') =2D t\ , 

k B Tn J 

with the short time diffusion coefficient Dq. Stress kernels generalize the (inverse) fluid 
isothermal compressibility n T to finite wave vectors and frequencies: G(q, t) = G(q, t) + 
(1/Kq )qq. Here, = S q /{kBTn) is given by the equilibrium fluid structure factor. The 
time-dependent stress kernels (with reduced dynamics according to the ZM-formalism 
[20] indicated by ig) 

N 

(25) G(q,t) = (n/fc B T)<a q (t e )V q > , with a q = -±= V F, 

reduce to the stress auto-correlation functions in the limit of vanishing wave vector, with 
(To a component of the potential part of the microscopic Irving-Kirwood stress tensor 
defined in Eq. (6): 

11 

(26) G(t) = Jim G q (i) = — (<7 (*)* ^o> • 

The shear element is measured as function of frequency lo in linear rheology, G'{lo) + 
iG"(uj) = iuj J °° dt e~ tu)t G(t), where G' is the storage and G" the loss modulus [4]. 

Considering glass a non-ergodic state, the stress kernels, like density (but not velocity) 
fluctuations [29], do not relax to equilibrium but take finite values at infinite time: 

(27) G(q,t -> ^^Goc (q), 

which shows up as finite zero- frequency elastic shear modulus in rheology, G'(ui — > 0) — > 
Goo. Equation (27) predicts using Eq. (24) 

(28) C (q,t^ & o)^C 00 (q)-2^(G 00 (q))- 1 . 

q 

Displacement differences stay below a finite limit for all times. Equation (22) is a valid 
approximation of displacement fluctuations in glass! The expansion e* q ' r< (*) = e tq-r< + 
C(q • Uj(t)) is justified, based on the collective variance as measure of the finite variance 
of individual particle displacements. This is an idealization as it requires that r Q >• At, 
which restricts the approach to states (far) away from the transition (deep) in the glass. 
Considering the local motion of the particles around the time-averaged positions, the 
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cquipartition theorem (3) follows and the dynamical matrix 2?(q) is given by the frozen- 
in stress fluctuations [28] 



Thus, all signatures of isotropic elasticity arise from treating glass as a metastable state 
with frozen-in stresses. The conservation law of momentum, which leads to the explicit 
factors q 2 in Eqs. (24,29), provides the route to explain the long-ranged nature of dis- 
placement fluctuations. 

Figure 4 shows 'dispersion relations' in a 2-dimensional colloidal system obtained 
by video microscopy [28]. Displacement fluctuations with time-averaged sites (see right 
panel), were correlated and the eigenvalues of the dynamic matrix in Eq. (3) diagonalized. 
Previous microscopy approaches had concentrated on the covariance matrix to obtain the 
density of states [31]. The dispersion relations display the signatures of a solid, viz. trans- 
verse elastic modes, in colloidal glass and in viscoelastic fluids at intermediate times; the 
latter aspect follows Maxwell's insight that fluids are solid-like at high frequencies [9]. 
Considering that stresses are notoriously difficult to measure in dense amorphous sys- 
tems, the smoothness of the dispersion relations obtained from the equipartition theorem 
is noteworthy. Computer simulations also explicitly tested relation (29), and determined 
the Lame coefficient \i and the frozen-in shear stress G m independently. Agreement, as 
predicted by Eq. (29) was observed within the statistical error bars [28]. 

5. Conclusions and outlook 

This short review compared different theories for the displacement correlations, elas- 
tic dispersion relations, and elastic moduli, which can be studied by techniques like video 
or confocal microscopy in colloidal solids. The effects of disorder, either as defects in pe- 
riodic solids, or in amorphous solids was highlighted. A number of theoretical relations, 
e.g. between fluctuation and DFT results, remain to be understood, and experiments 
in three-dimensional and two-dimensional solids would be useful. Extensions to optical 
phonons in periodic solids, and phasonic modes in quasi-crystalline ones [32] are promis- 
ing. In glasses, the role of defect motion remains open, and extensions to capture the 
Boson peak and the density of states [33] should be included in the theoretical approach. 

Appendix A. 

The free energy density / = F/V is derived starting from the free energy relation 
in Eq. (1) in order to determine its dependence on the thermodynamic fields; one finds 
/ = fi u afi,n), respectively / = f(u a p,c), with u a p the strain, n = N/V the total 
density, and c the defect density [3, 2]. Constant temperature (dT = 0) is considered. 

Because of the homogeneity of F in the thermodynamic limit, one obtains 



(29) 




/ = F/V = [in + hap u a p - p(n, h a p) . 
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Here the Gibbs-Duhem relation was already indicated, which states the dependence of 
the pressure on chemical potential p and stress tensor h: 

dp = n dp, + u a p dh a p . 

It follows from F — pN — Tr[h-u] = — pV being extensive with V as single extensive 
variable. Combining the total differential of the first equation with the second gives 

df = pdn + h a[i du a j3 , 

the dependence of the free energy density on strain and total density. Transforming to 
defect density c via dn = —ndu aa — dc, gives the free energy density as function of strain 
and defect density: 

df = —p dc + a a p du a [j , with a a p = h a p — np S a p . 

Fleming and Cohen [22] rewrite / in terms of the free energy density / = F/N per 
particle using df = d [ nfj = ndf + f dn. This leads to 

P 

ndf = —- dc + t a f! du a/3 , with t a p = h af} - (np - f) 5 a p . 

The stress tensor t gives the reversible part of the momentum conservation equation 
[22, 7], which shows that it can be obtained from averaging the microscopic Irving- 
Kirkwood stress tensor in Eq. (6). For fixed N, the last equation recovers Eq. (5). In 
unstrained equilibrium, this confirms Eq. (7). 

* * * 
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